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INTRODUCTION 

The results to be reviewed are divided into two categories : those that ref- 
late two-dimensional shell theory to three-dimensional elasticity theory and 
those concerned with shell theory per se . In the second category I further dis- 
tinquish between results for general elastic systems that carry over, by spe- 
cialization or analogy, to shells and results that are unique to shell theory 
itself. Because of the limitations of space and my interests, I do not men- 
tion multilayered or sandwich shells. A good discussion of these with an 
ample list of references may be found in Librescu's book [l]. Also, in view of 
the excellent review articles by Stein [2] and Hutchinson and Koiter [S], I 
have not attempted to review the enormous literature on shell buckling. 

TWO APPROACHES TO SHELL THEORY 

Most texts derive shell theory by a mixture of two-and three-dimensional 
considerations. However, a number of recent papers have adopted one of the 
following two extreme approaches: 

A. A shell is idealized as a material surface in three-dimensional Eu- 
clidean space capable of transmitting forces and moments. The physical laws 
for this two-dimensional eontinutim are postulated in analogy with those for a 
three-dimensional one. Stress-strain laws and even failure criteria are for- 
mulated in terms of two-dimensional variables and may be deduced directly from 
experiments on the shell material. The papers by Sanders [i|], ’Ericksen and 
Truesdell [5], Serbin [6], Budiansky [ 7 ], Simmonds and Danielson [8], and 
Reissner [ 9 ], to mention but a few, as well as much of the monumental treatise 
by Naghdi [lO] are written in this spirit. 

B. No matter how thin, a shell must be regarded as a three-dimensional con- 
tinuum. However , the governing equations can be enormously simplified by con- 
sidering various formal asymptotic expansions of the unknowns in terms of ap- 
propiate "thinness’’ parameters. In the interior of the shell (i.e. away from 
edges, concentrated loads or geometric discontinuities of one sort or another) 
the leading terms of the expansions satisfy various sets of two-dimensional 
equations that we call, collectively, the shell equations. Among those who 

have contributed recently to this second approach are Green [ll], Johnson and 
Reissner [12], Cicala [13], Van der Heijden [l^], and especially Goldenveiser 
(see the references cited in [15].) 
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The virtue of the first approach is also its shortcoming: there is no way 
to estimate intrinsically the errors made hy negleeting three-dimensional ef- 
fects. Or, from another viewpoint, there is no systematic way to construct a 
refined shell theory. 

A drawback of the second approach, aside from its tediousness, is that it 
requires a knowledge at the edges of the shell of the distribution in the 
thickness direction of the applied stresses or displacements, As Koiter has 
emphasized [l6], we never know these distributions precisely, except at a free 
edge. Another drawback of the second approach is that, because the thickness 
of the shell is always incorporated in the expansion parameters, one set of 
uniformly valid interior (i.e. shell) equations does not emerge. Rather there 
is one set of equations for a ''membrane" state-, another for an "inextensional 
bending" state, another for a "simple edge effect", another for a "degenerate 
edge effect", and, if one is dealing, for example, with an infinite cylindrical 
shell subject to self-equilibrating edge loads, still another set of equations 
is needed to recover the "semi -membrane" theory of Vlasov [l7j P* 25^+]. 

THE ASYMPTOTIC APPROACH 

The goal here is to provide a systematic method of refining the analysis 
of thin-walled bodies. One important consequence of the asymptotic approach 
is the verification and refinement of the classical Kirchhoff boundary condi- 
tions. Another useful result is that it gives a method for computing the do- 
minant stresses in the immediate vicinity of an edge without the need of 
solving a full three-dimensional problem. We shall first illustrate the es- 
sence of . the asymptotic method by means of a simple example drawn from the 
work of Goldenveiser and Van der Heijden. Then we shall indicate the implica- 
tion of the results for nonlinear shell theory. 

Let (r,9,z) denote a set of cylindrical coordinates and consider a homo- 
geneous, elastically isotropic plate that occupies the region O^r^R, -H<z^H. 
Let the plate be free of body forces and edge tractions but subject to self- 
equilibrating normal tractions on its upper and lower faces. The linear 
equations of elasticity may be expressed as three equilibrium equations for 
the six Independent components (c > T, cfg,t ,Tg,a) of the symmetric stress 
tensor plus six stress-strain relations witK the strains expressed in terms of 
the components (u,v,w) of the displacement vector. Let p=r/R and ?=z/H. Then 
the boundary conditions read 

a( P , 0 , + 1 )=+ ( P , e ) , ( P , 0 ,+1)=7 q ( p , 0 ,± i )=0 ( 3. i ; 

a^(l,0,5)=T(l,0,;)=t^(l,0,C)=O, (3.2) 

where is a reference stress chosen so that |p|^ 1. The boundary conditions 

induce a state of pure bending in which (a^jT^ag jCr^u^v) are odd in C and 

(t jT^jw) are even, 
r 6’ 

Goldenveiser Vs approach, following earlier work by Friedrichs and Dressier 
[l8] and Green [ll], is to express each unknown as the sum of a "basic” or in- 
terior contribution plus two distinct "auxiliary" or e^ge zone contributions. 
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The edge zone contrihutions are expressed in terms of the sealed yariahle 
?=(R-r)/H=e"^(l-p) SO that 5 for example, 

( P 9 9 9 C 9 ^ ) ” ^Q t ^ P 5 ® 9 C 9 ^ ( C 9 9 5 C 9 ^ C 9 9 5 C ^ £ ) ] 

and ■u(p,e,?;e)+(B/E)[U(p, 0 ,C;e)+u(C,e, 5 ;e)+u(^, 0 ,C;e)]. 

For a traction free edge, Goldenveiser [I9] assumes the following formal as- 
ymptotic expansions 


(Z^,T,Eq,T^,Tq,E,U,V,¥) 


“fe'' ( ,t“e“ ,eT“,eT“ e^E“ ,u“ ,v“ ,e" V ) 


- 1 , 


(er^, > . . ,w;~2e (€a^,x ,eaQ,eT^,Tg,e0 ,e u,ev,e w) 

(8, «)-e|e“(8j,rt",8^ 

The edge zone contributions are assumed to vanish exponentially as 


(3.3) 

( 3 .i^) 

( 3 . 5 ) 


When these representations are suhstituted into the elasticity equatiohs 
and their assimied asymptotic character accounted for, there results an infi- 
nite sequence of differential equations for each infinite sequence of coeffi- 
cients . . . ,W^},{o^». . . Furthermore, the boundary con- 

ditions (3.1) and (3.2) imply that for C=+d, 


and that for p=l and C= 0 ; 

(j“,T°+T°y+T^,£"'"^+3®+8“ T“'"2+T”'"^+t”,T”+T”+t”)=0 

r’ ’ ’ r r r, ’ r r r 


(3.6) 

(3.7) 

(3.8) 


where n= 0 ,l, 2 ,,.. . 

The equations for the interior coefficients may be integrated systemati- 
cally with respect to Application of the face boundary conditions ( 3 . 6 ) 
leads, in the first instance, to the classical equation of plate bending 


( 2/3 )MW°= (1-V^ )p , AW°=¥°^^+W*^^^ . 


( 3 . 9 ) 


All of the remaining lowest order interior coefficients are expressible in 
terms of ¥° ; in particular 


E°=-(l-V^ )"^C [¥°pp+v( p“^ 7 °p+p"^°gg JE 2 / 3 CM® (p , ) 


( 3 . 10 ) 
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(3.11) 

( 3 . 12 ) 


T°= (1+v [ p"^°Q-p“ Vp 0 ] =2/3CH° ( p , e ) 

T°=-|(l-V^)"^(l-C^)(AW°),p= 3/4 (l-C^)Q°(p,e). 

The first of the edge boundary conditions in (3.8) , namely yields only 

one of the tvo boundary conditions needed for ¥ . To obtain the second, the 
edge zone solutions must be considered. 

The infinite sequence of differential equations for the set of edge zone 
coefficients (a can be grouped into sets i^Thich resemble the nonho- 

mogeneous St.-Venant equations for the torsion of a prism whose cross-section 
is the semi-infinite strip Likewise the differential equations for 

the coefficients w) can be grouped into sets which resemble the non- 

homogeneous equations of plane strain for the same semi-infinite strip. The 
solutions of the torsion and plane strain problems are coupled through the non- 
homogeneous terms in the differential equations as well as through the boundary 
conditions ( 3 . 7 ) and ( 3 . 8 ) which also link these solutions with the interior 
solutions. It should be noted that in the edge zone differential equations, 6 
appears only as a parameter. 

In order that the edge zone solutions decay as it is necessary that 

the forces and moments applied to the boundary of the semi-infinite strip be 
equilibrated by the non-homogeneous terms in the torsion and plane strain 
equilibrium equations. These integral conditions yield, ultimately, the addi- 
tional boundary conditions needed for the various interior solutions. For ex- 
ample, the Kirchhoff boundary condition that relates the shear stress resul- 
tant and the derivative along the edge of the twisting stress couple H is 
obtained as follows. 


The solution of the lowest order torsion problem may be expressed in terms 


of a stress function ijj , where 


Aip°= 0 , 1 |J,^( 0 ,C)=-C 

(3.13) 

and 


T°= -2/3CH°(1,0)t{)°^, t°= -2/3CH°(l,e)il;°^ 

(3.14) 


The lowest order equation for equilibrium in the ^-direction for the plane 
strain problem is 




(3.15) 


From the last of the hoxindary conditions (3.T) and ( 3 . 8 ), the condition that 
the net forces ia the ^-direction add to zero, to lowest order, is 

[T°(i,e,c) + /” R°(c,e,c)de]dc=o. (3.16) 

With the aid of (3.12) and (3.13) to (3.15), (3.l6) reduces to 

Q° + H° = 0 at p=l, (3.17) 

r u 


620 



which is the second boundary condition for W°. It is important to note that 
one never needs to actually solve for Jp to obtain this result. 


GOLDENVEISER’S EXTENSION AND KOITER’S SIMPLIFICATION OF THE PRECEDING RESULTS 

The solution for (E^, ... , reduces to the solution of a biharmonic 

equation for To obtain boundary conditions for one again considers the 

integral conditions of overall equilibrium necessary to guarantee decaying 
edge zone solutions. To evaluate these, one must solve explicitly for ^ 
(which is easily done) but needs only to consider the form of the solution of 
the lowest order plane strain problenr. After a straightforward but tedious 
analysis, there results the refined boundary conditions of • Goldenveiser [19]: 


mJ = AH,°, + AH°g=0, (4.1) 

where A=1.260,.,is computed from the solution for The details of the cal- 

culations leading to (4.1) may be found in a report by Van der Heijden [20], 


Goldenveiser ’s results may be restated in the following useful way. Con- 
sider a plate of radius R and thickness 2H subject to a self-equilibrated nor- 
mal pressure p but otherwise free of surface and edge tractions. Solve the 
classical equation of plate bending subject to the refined boundary conditions. 

- A(H/R)H,g = 0, + H,g + A(H/R)H,g = 0 at p=l. (4.2) 

Then the stresses in the interior of the plate, to within a relative error of 
0 (r2/r 2), are given by the formulas for T°, etc. but with W replaced by 
W. Moreover, in the edge zone of the pla^e, the dominant stresses, to within a 
relative error of 0(H/R), are given by these same formulas except that T is 
replaced by T°+x , and Tg is replaced by Xg, where x and Xg are given by (3.l4) 

These results are simple and satisfying. Though derived for, perhaps, the 
simplest, non-trivial problem imaginable, their qualitative implications for 
shells with free edges undergoing large deformations is clear, namely l), the 
most importnat refinement of the classical shell equations are in the boundary 
conditions and 2), the dominant stresses near a free edge can be inferred from 
the solution of the shell equations and the solution of a torsion problem for 
a semi-infinite strip. To give these statements a quantitative form via an 
asymptotic analysis would seem to be a formidable task. 


The problem of refining the Kirchhoff boundary conditions at a free edge 
has, fortunately^, been solved by Koiter [15] in an alternate way, using an in- 
genious energy argument. As Danielson [21], and Koiter [22] have shown, the 
three-dimensional tangential shear stress predicted by shell theory -at a free 
edge does not vanish, even though the Kirchhoff botindary conditions are satis- 
fied exactly. Thus the conventional strain energy expression of shell theory 
overest imates the torsional energy in the neighborhood of a free edge. To as- 
sess this error, Koiter considers the torsional rigidity of a flat strip whose 
thickness is equal to that of the shell. By comparing this expression with 
that given by classical plate theory he is able to identify an edge zone cor- 
rection factor which is proportional to the twist per unit length of the edge 
of the strip. The torsional energy associated with this term is therefore 
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expressible as a line integral. For an arbitrary shell with a smooth edge 
curve Koiter argues that one merely needs to insert an appropiate expression 
for the edge twisting per unit length for the shell into this line integral and 
then subtract this expression from the conventional surface integral for the 
shell energy. 

Koiter' s result may be of limited practical value. If the shell has other 
edges that are not free of stress, it is most likely that the associated shell 
boundary conditions cannot be refined because the corresponding boundary con- 
ditions of elasticity theory cannot be determined precisely. The shell eq^ua- 
tions are elliptic, hence the influence of boundary conditions extend every- 
where, and it would be inconsistent to use refined boundary conditions at a 
free edge but unrefined ones at another edge. 

The results of this section also imply that so-called thick shell theories 
are meaningless if applied to homogeneous shells with edges. ¥e should note, 
however, that Van der Heijden has shown that Reissner's latest thick plate 
theory [23] does give fairly good nmerical results for stress concentration 
factors for circular holes in infinite plates. 


THE DIRECT APPROACH TO SHELL THEORY' 


Here and in the following section we mention briefly — space limitations 
permit no more ■— some recent work concerning different formulations, implica- 
tions, simplifications and the reduction of certain problems of the now gen- 
erally accepted equations of first-approximation shell theory. 


Formulations of the Nonlinear Theory 

A strictly mechanical theory of shells may be expressed entirely in terms 
of the midsurface displacement components [5 ]. If dynamic effects are occluded ,. 
alternate formulations are possible in terms of the components- of a stress 
function and rotation vector [8], or in terms of stress resultants and bending 
strains [l5l. In the last case, any displacement boundary conditions need to 
be reformulated in terns of strains [24,25]. This in itself has advantages, 
for it automatically leads to the boundary conditions for inextensional defor- 
mation and, in the linear theory, it gives boundary conditions that are the 
geometric analogues of the Kirchhoff conditions. 

Thermodynamic Considerations 

These are important for at least three reasons, l) heating a shell may 
cause it to fail, buckle, or vibrate; 2) the best justification of the static 
approach to stability for a continuous body is a thermodynamic one; and 3) the- 
coupling of mechanical and thermal effects produces damping. 

There is a plethora of papers on l) that we shall not attempt to review; 
a few texts give a discussion of the underlying ideas. The thermodynamic as- 
pects of stability in general elastic systems are discussed in [ 26 , 27 , 28 ]. 

These results are direetly transferable to shell theory. The specific form 
and role of the laws of thermodynamics in shell theory are discussed. in [lO]. 
The effect of thermal damping on the free vibrations of shells is considered 
in [ 29 ] where it is also shown that, because the damping is light, perturbation 
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methods may he used to advantage. 

Variational Principles 

A problem of long standing in nonlinear elasticity has been to. formulate 
a principle of complementary energy. Recent work [30,31,32] has established 
conditions under which this is possible. In particular, in [33] and [3^], 
these resiults have been applied to the nonlinear von Karman plate eq.tmtions 
and Marguerre shallow shell equations to obtain upper and lower bounds on an 
associated energy fimctional, 

SOME HEW RESULTS IH LIHEAR SHELL THEORY 
Shells As Beams 

For general cylindrical shells and shells of revolution, one may consider 
special classes of solutions that, in a St. Venant sense, , correspond to the 
stretching, bending, twisting, and flexure of a beam. In many cases the re- 
sulting equations can be solved explicitly. See [35j36], 

Reduction of the Governing Equations 

The shell equations constitute a system of eighth order. For analytical 
purposes, especially for the application of pertiirbation methods, it is often 
convenient to attempt to express these eqxoations as two coupled fourth order 
equations. (A single eighth order equation destroys the very useful static- 
geometric duality). Such reductions have been found for spherical, general 
cylindrical, and minimal shells as well as for shells of revolution. A reduc- 
tion for arbitrary, non-developable shells is also possible, but does involve 
some loss of accuracy. See [37l where other references are cited. 

Membrane Theory 

It is well known that shells with the proper shape and boundary support 
can be analyzed with good accuracy by monbrane theory. The details of such 
an approach are spelled out in a very general but useful way in [38], 

Cracks and Cutouts 

Shells may contain cutouts by design and cracks by accident. In practice 
the dimensions of these cracks and cutouts is apt to be small conq)ared to some 
characteristic geometric dimension of the shell, permitting shallow shell the - 
ory to be applied. The calculation of the stresses has been reduced to the 
solution of coupled singular integral equations [39] that have been solved 
numerically for several important problems. See [J+O] and the references cited 
there . 

Pointwise Estimates For Approximate Solutions 

The Prager-Synge hypercircle method is useful for constructing approxi- 
mate solutions to linear- shell problems, and provides mean square error es- 
timates for the approximate stress field. More desirable are pointwise es- 
timates for both the approximate stress field and the approximate displace- 
ment field. For recent work on this problem see [i^l] and the references cited 
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therein. 


¥ave Propagation, Asymptotics, and St-Venant’s Principle 

These are three additional areas in vhich there has been significant re- 
cent progress but which cannot be reviewed for lack of space. 
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